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The simplest, "standard" model of Big Bang Nucleosynthesis (SBBN) assumes three light neu- 
trinos (N„ = 3) and no significant electron neutrino asymmetry (y e — P e asymmetry parameter 
£ e = n e /kT, where /x e is the v e chemical potential) leaving only one adjustable parameter: the 
baryon to photon ratio n = ns/n-f. The primordial abundance of any one nuclide can, therefore, 
be used to measure r\ and the value derived from the observationally inferred primordial abundance 
of deuterium closely matches that from current non-BBN data, primarily from the WMAP survey. 
However, using this same estimate, there is a tension between the SBBN-predicted helium-4 and 
lithium-7 abundances and their current, observationally inferred primordial abundances, suggesting 
that N„ may differ from the standard model value of three and/or that £ e may differ from zero (or, 
that systematic errors in the abundance determinations have been underestimated or overlooked). 
The differences are not large and the allowed ranges of the BBN parameters (r), N v , and £ e ) per- 
mitted by the data are quite small. Within these ranges, the BBN-predicted abundances of D, 3 He, 
4 He, and 7 Li are very smooth, monotonic functions of 7710, AN V = N„ — 3, and £ e . As a result, it 
is possible to describe the dependencies of these abundances (or powers of them) upon the three 
parameters by simple, linear fits which, over their ranges of applicability, are accurate to a few 
percent or even better. The fits presented here have not been maximized for their accuracy but, 
rather, for their simplicity. To identify the ranges of applicability and relative accuracies, they are 
compared to detailed BBN calculations; their utility is illustrated with several examples. Given the 
tension within BBN, these fits should prove useful in facilitating studies of the viability of various 
options for non-standard physics and cosmology, prior to undertaking detailed BBN calculations. 



I. 



BIG BANG NUCLEOSYNTHESIS 



Shortly after the emergence of the Big Bang model it 
was realized that conditions were ripe during the early 
Universe for a brief period of nucleosynthesis. Just as 
in every other setting where nuclear reactions occur, the 
yields of the emerging nuclides are governed by three en- 
vironmental characteristics: the duration of the event, 
the density of the reactants, and their thermal proper- 
ties. A brief, dilute and cool environment would yield a 
very different set of abundances compared to those that 
would emerge if Big Bang Nucleosynthesis (BBN) had 
been long-lasting, dense and hot. Although it is possible 
to characterize BBN in such broad generality (with some 
success; see, for example, [1, 2]), the paradigm most fre- 
quently encountered uses the Friedman equation to relate 
the BBN expansion rate, H, to the thermal properties of 
the particles present at that epoch 
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with Gn being Newton's constant and p the total energy 
density. The standard model of particle physics provides 
the candidate particles whose energy density contribute 
to p, but it is possible that this may fall short of the actual 
particle content at that time. Even so, after blueshift- 
ing the currently observed Cosmic Background Radia- 
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tion (CBR) energy density back to the epoch of BBN, 
it emerges that the energy density of the Universe was 
dominated during BBN by the relativistic particles. The 
thermal properties of each particle species are described 
by a temperature and a chemical potential, but thermal 
coupling among the constituents equilibrates the temper- 
atures while chemical coupling relates the chemical po- 
tentials. Charge neutrality provides additional relations. 
This reduces the number of independent quantities with 
the result that BBN can be described by a minimal set 
of parameters, typically taken to be: the density of the 
nucleons/baryons, the chemical potentials of the three, 
active neutrinos, and the energy density in unaccounted 
particles (and/or additional terms on the right hand side 
of eq. 1). In what follows the role of each is discussed and 
the parameters used to describe them are introduced. 



A. The Baryon Density: r\ 

The simplest, "standard" model for BBN (SBBN) sets 
any "extra" energy density and the chemical potentials of 
the three neutrinos to zero leaving just the density of the 
baryons, jib, as the only adjustable/free parameter. But 
since ne drops as the Universe expands, it is useful to 
introduce the baryon to photon ratio r\ = n^/n^, (7710 = 
10 10 ?7), since this quantity is dimensionless and constant, 
except during the period of e annihilation. 

From inspection of the flow of nuclei through the reac- 
tion network in SBBN it can be seen that the primordial 
abundances of the relic light nuclides D, 3 He, and 7 Li 



are determined by the competition between the nuclear 
production and destruction rates, which scale with the 
nucleon density. To an extent which varies from nuclide 
to nuclide, these three nuclei are all potential baryome- 
tcrs. In contrast, as the light nucleus with the largest 
binding energy per nucleon, the 4 He relic abundance is 
relatively insensitive to the magnitudes of the nuclear re- 
action rates because, for the range in baryon density of 
interest here, they are sufficiently rapid to burn nearly 
all neutrons present initially to 4 He. What sensitivity it 
does possess is due to the role that r\ plays in determin- 
ing the temperature at which the D abundance becomes 
significant: it is the build up of a large D abundance 
that leads to the most rapid phase of nuclear burning 
in BBN and the termination of free neutron decay as an 
important process. 



B. Early-Universe Expansion Rate: S 

The early Universe is radiation-dominated, with its ex- 
pansion rate determined by the energy density in the 
relativistic particles, pr. Prior to BBN, e.g., at a tem- 
perature of a few MeV, before e annihilation, the stan- 
dard model of particle physics provides photons, e pairs 
and three flavors of left-handed (i.e., one helicity state) 
neutrinos (and their right-handed, antineutrinos) as con- 
stituents of this dominant component. With all chemical 
potentials set to zero the energy densities are related by 
thermal equilibrium so that 



PR = Pi 
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where p 7 is the energy density in the CBR photons 
(which, today, have redshiftcd to become the cosmic 
background radiation photons observed at a temperature 
of 2.7K). In this case, the time-temperature relation (de- 
rived from the Friedman equation) is, 

Pre - e ± annihilation : t T^ = 0.738 MeV 2 s. (3) 

In SBBN it is often assumed that the neutrinos are 
fully decoupled prior to e annihilation and that they 
don't share in the energy transferred from the annihilat- 
ing e pairs to the CBR photons. In this approximation, 
the photons are hotter than the neutrinos in the post- 
e ± annihilation universe by a factor T^/T u = (11/4) 1 / 3 , 
and the relativistic energy density is 

PR = Pi + 'iPv = l-68p 7 , (4) 

corresponding to a new time-temperature relation, 

Post - e ± annihilation : t T* = 1.32 MeV 2 s. (5) 

In quite general terms, one possible consequence of 
new physics beyond the standard model could be a non- 
standard, early Universe expansion rate, parameterized 
by an expansion rate factor S, 



Although the introduction of S here is quite general, in 
practice one must adopt a specific scheme through which 
this occurs. Such a non-standard expansion rate might, 
but need not, be due to a change in the energy density: 
a change in the strength of gravity would also alter the 
expansion rate of the early Universe [3] as would non- 
standard, higher dimensional models which modify the 
expansion rate - energy density relation (the Friedman 
equation, eq. 1) in our 3+1 dimensional world [12]. Dif- 
ferent gravitational couplings for fermions and bosons [4] 
would have similar effects. These different mechanisms 
for implementing a non-standard expansion rate are not 
necessarily equivalent. 

If consideration is restricted to only the possibility of 
additional energy density, then 
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where p' K = pr + px and X identifies the extra com- 
ponent. With the further restriction that the X arc 
also relativistic, this extra component behaves just like 
an additional neutrino though we emphasize that "X" 
need not be additional flavors of active or sterile neu- 
trinos. In these circumstances S is a constant prior to 
e annihilation and it is convenient to account for the 
extra contribution to the standard-model energy density 
by normalizing it to that of an "equivalent" neutrino fla- 
vor [11], so that 



7 
px = AN vPv = -AN vPl 



For this case, 
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However the expansion rate is more fundamental to BBN 
than is A7V„, so we parameterize this class of non- 
standard models using S (but, for comparison, we will 
often also quote the corresponding value of AN V as given 
by eq. 9). 

Not every additional energy density component can be 
accommodated this way since we are requiring that px is 
proportional to p v as the Universe evolves. An example 
that breaks this mold is if the "X" were not 'decoupled', 
in the sense that X shared in the energy released when 
the e pairs annihilate. Other examples are found within 
some quintessence models [3, 5, 6, 7]. 

As we discussed above, D, 3 He, and 7 Li act as the 
principle baryometers for BBN since the 4 He relic abun- 
dance does not vary considerably with n. But 4 He is very 
sensitive to the competition between the weak interac- 
tion rates (interconverting neutrons and protons) and the 
universal expansion rate which, during the early, radia- 
tion dominated evolution is fixed by the energy density 
in relativistic particles ("radiation"). As a result, while 
D, 3 He, and Li probe the baryon density, consistency 
between their abundances and that of 4 He tests the stan- 
dard model and provides constraints on (and, perhaps, 
hints of) physics beyond the standard model. 



C. Lepton Asymmetry: £ 

The baryon-to-photon ratio rj = n-g/n^ provides a di- 
mensionless measure of the universal baryon asymmetry 
which is very small (rj ^ 10 -9 ). By charge neutrality the 
asymmetry in the charged leptons must also be of this 
order. However, there are no direct observational con- 
straints (see [10, 13, 14] and further references therein), 
on the magnitude of any asymmetry among the neutral 
leptons (neutrinos). A dimensionless measure of the mag- 
nitude of the neutral lepton asymmetry is provided by £, 
the ratio of the neutral lepton chemical potential to the 
temperature (in energy units): £ = fi/kT. For exam- 
ple, for a neutrino flavor a, the asymmetry ("neutrino 
degeneracy" ) L a , between the numbers of v a and v a is 



L a = 



12C(3) \?* 7T 2 
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Note that for £ Q <C 1, L a fa 0.684£ Q . Mixing among 
the three active neutrinos (z/ e , v^, v T ) ensures that at the 
time of BBN, L e « L M w L T (& » £„ » &.) [15]. 

While any neutrino degeneracy (£ Q < as well as > 0) 
will increase the energy density in the relativistic neu- 
trinos and, hence, the early Universe expansion rate, the 
discussion here is limited to sufficiently small values of £ Q 
so that the effect on S of a non-zero £ Q is negligible. Even 
so, a non-zero but relatively small asymmetry between 
electron type neutrinos and antineutrinos (£ e ^ 10 -2 ), 
while large compared to the baryon asymmetry, can still 
have a significant impact on the early Universe. The in- 
terconversion of neutrons to protons is affected by the v e 
so that any non-zero £ e will alter the n/p ratio, thereby 
affecting the yields of the light nuclides formed during 
BBN. 

Of the light nuclei, the neutron limited 4 He abundance 
is the most sensitive to £ e and, together with the D, 3 He, 
and 7 Li baryometers, again provides a test of the consis- 
tency of the standard model. 



D. BBN Predictions 

The BBN-predicted D, 4 He, and 7 Li isoabundance 
curves for S versus rj and for £ e versus n are shown in 
Figures 1 and 2. For D, the curves represent the ratio (to 
hydrogen) by number y-Q = 10 5 (D/H) = 2, 3, 4; similarly, 
for 7 Li, the curves correspond to yu = 10 10 (Li/H) = 3, 
4, 5; for 4 He, the curves correspond to the mass fraction 
Y P = 0.23, 0.24, 0.25. These figures show clearly that D 
and Li are good baryometers, sensitive to 77, while He 
is a good chronometer, sensitive to S, as well as a good 
leptometer, sensitive to £ e . 



II. OBSERVED PRIMORDIAL ABUNDANCES 

The success of Big Bang Nucleosynthesis relies on its 
ability to predict the observed primordial abundances 




FIG. 1: BBN-predicted isoabundance curves for D (dotted), 
7 Li (dashed), and 4 He (solid) in the expansion-rate parameter 
(S), baryon density parameter (77) plane. From left to right, 
for D, y D — 4, 3, 2, while for 7 Li, y L i = 3, 4, 5. For 4 He, from 
bottom to top, Y P = 0.23, 0.24, 0.25. 
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FIG. 2: BBN-predicted isoabundance curves for D (dotted), 
7 Li (dashed), and 4 He (solid) in the lepton asymmetry pa- 
rameter (£ e ), baryon density parameter (rj) plane. From left 
to right, for D, j/d = 4, 3, 2, and for 7 Li, j/Li = 3, 4, 5. For 
4 He, from top to bottom, Y P = 0.23, 0.24, 0.25. 



and, conversely, to learn about the cosmological param- 
eters using these same abundances. However the com- 
parison of predictions and observations is far from trivial 
because of the further processing of the nuclei since the 
end of BBN. 

Deuterium is the baryometer of choice, primarily be- 
cause since BBN its observed abundance should have 
only decreased as gas is cycled through stars where D 
is burned to heavier, more tightly bound, nuclei [18]. As 
a result, D observed anywhere in the Universe, at any 
time during its evolution, should provide a lower bound 
to its primordial abundance. Further, the deuterium ob- 
served in the high redshift, low metallicity QSO absorp- 



tion line systems (QSOALS) should be very nearly pri- 
mordial. In contrast, the post-BBN evolution of 3 He and 
of 7 Li are considerably more complicated, involving the 
competition between production, destruction, and sur- 
vival. As a result, at least so far, the current, locally ob- 
served (in the Galaxy) abundances of these nuclides have 
been of less value than has that of deuterium. Indeed, 
over the range in j/d to be adopted below, the primor- 
dial abundance of 3 Hc is predicted to lie in the narrow 
range, 1.0 < 7/3 < 1.1, in excellent agreement with that 
inferred from Galactic observations [19]. Thus, 3 He pro- 
vides similar, but less compelling constraints than does 
D. While current estimates of primordial lithium, based 
on observations of very metal-poor (nearly primordial) 
stars, arc problematic, 7 Li is retained here in order to 
highlight this challenge. For 4 He the post-BBN evolu- 
tion is straightforward, but difficult to quantify at a suf- 
ficient level of accuracy. As gas is cycled through stars, 
hydrogen is burned to helium; only a very small fraction 
of the helium is burned to heavier nuclei. If material is 
sampled "here and now" , a correction based on models 
of the chemical evolution of galaxies would need to be 
made for post-BBN produced 4 He. To avoid the model- 
dependent uncertainties associated with such corrections, 
which would surely overwhelm the observational errors, 
the best strategy is to search for 4 He in the least evolved 
objects, in this case the metal-poor, extragalactic H 11 re- 
gions. 

Inferring the primordial D abundance from the 
QSOALS has not been without its difficulties, with some 
of the original abundance claims having been withdrawn 
or revised. Presently there are 5-6 QSOALS with rea- 
sonably firm deuterium detections [20, 21, 22, 23, 24, 25]. 
However, there is significant dispersion among the abun- 
dances and the data fail to reveal the anticipated "deu- 
terium plateau" at low metallicity or at high redshift [26] . 
Here we adopt the five abundance determinations col- 
lected in the recent paper of Kirkman et al. [25]. The 
weighted mean value of j/d is 2.6 . The dispersion 
among these five data points is very large, resulting in 
a x 2 = 15.3 for four degrees of freedom, suggesting that 
one or more of these abundance determinations may be in 
error, or affected by unidentified and unaccounted for sys- 
tematic errors. To allow for this, we follow the approach 
advocated by [22] and [25] and adopt for the uncertainty 
in j/d the dispersion among the data points divided by 
the square root of the number of them. Thus, the primor- 
dial abundance of deuterium to be used here is chosen to 
be: j/d =2.6 ±0.4. 

A somewhat less than clear situation also exists for 
the determinations of the primordial abundance of 4 He. 
There have been two, largely independent, estimates of 



1 This differs from the result quoted in Kirkman et al. because 
they have taken the mean of log(j/o) and then used it to infer 
j/ D (y D = 10< ; °^d)>). 



Yp based on analyses of large data sets of low- metallicity, 
extragalactic Hn regions. The "IT" [27, 28] estimate 
of Y(IT) = 0.244 ± 0.002, and the "OS" determination 
[29, 30, 31] of Y(OS) = 0.234 ± 0.003 which differ by 
nearly 3cr. Very recently, IT have expanded their original 
data set and attempted to account for some (but not 
all) of the systematic uncertainties [32]; see, e.g., [33]. 
For their full data set of 89 H11 regions, IT find Yp = 
0.2429 ± 0.0009. To provide a contrast with previous 
results, this value is adopted for our analysis here. 

As with 4 He, the abundance of 7 Li grows in the course 
of post-BBN evolution, with lithium being produced in 
some stars (at some distinct times in their evolution) as 
well as in cosmic ray spallation reactions (where CNO nu- 
clei are broken down to, among other nuclides, lithium). 
So, as for D and 4 He, lithium should be observed in the 
least evolved, most metal-poor objects. While lithium 
has been observed in the Sun and in the solar system, as 
well as in the interstellar medium of the Galaxy, these 
provide evolved, metal-enriched samples of little use in 
estimating the primordial abundance of lithium. The 
most relevant data comes from studies of the very metal- 
poor stars in the halo of the Galaxy or in Globular Clus- 
ters. While this material should be very nearly primor- 
dial, it must be kept in mind that these are observations 
of the surfaces of the oldest stars in the Galaxy. If, dur- 
ing the course of their evolution, the surface material of 
these stars were mixed with the lithium-depleted interi- 
ors, the currently observed surface lithium abundances 
may not reflect the lithium abundances in these stars at 
their birth. The most recent data from studies of halo 
stars is from Ryan et al. [34] who find yu = l-23lg' 3 6 . In 
contrast, and in some conflict with this result, Bonifacio 
et al. [35] derive from a sample of Globular Cluster stars, 

,</Li — z - la -0.38- 

Below we compare each of these estimates of the pri- 
mordial abundances with the BBN predictions, standard 
as well as non-standard, fixed by the adopted D and 4 He 
abundances. 



III. FITS 

The non-linear, coupled differential equations of BBN 
are not conducive to analytic solution so that detailed 
comparison of the theoretical predictions with the ob- 
served/derived abundances of the light nuclei can only 
be achieved after resorting to a lengthy process of nu- 
merical calculations. This necessity blurs the connection 
between the parameter set {77, S, £ e } and the data set 
{yu, Yp, yLi}, and, conversely, the parameter values the 
data recommend. Yet it is clear from Figures 1 and 2 that 
the relic light nuclide abundances are smoothly varying, 
monotonic functions of 77, S, and £ e , so it isn't surpris- 
ing that over limited but substantial ranges, simple rela- 
tions exist between the predicted abundances and these 
parameters. While the BBN-predicted primordial abun- 
dances are certainly not linearly related to the baryon 



density, the expansion rate, or the lepton asymmetry, 
our goal here is to find simple, linear fits to the predicted 
abundances (or powers of them) as functions of these pa- 
rameters (to be defined more carefully below). These 
fits work well, sometimes remarkably well, over limited 
ranges in the parameters (and/or over limited ranges in 
the predicted abundances of D, 4 He, and 7 Li). 

Since the domains over which our fits are applicable 
are restricted, we must focus upon specific values for 77, 
S and £ e . The "target" value/range of the baryon density 
is motivated by the (non-BBN) results from the WMAP 
constraints from the CBR temperature fluctuations [8] 
where, for SBBN, 7710 = 6.14 ± 0.25. As shown in [9, 10] 
(and references therein), this estimate is little affected by 
a non-standard expansion rate and/or any (small) lepton 
asymmetry. As a result, a limited range in 77 10 , centered 
around 7710 = 6, is chosen here: 4 < 7710 < 8. For our 
simple, linear fits to the BBN-predicted abundances as 
a function of S to be sufficiently accurate, we must re- 
strict the range in S. This is the case for the fits we 
adopt provided that 0.85 < S < 1.15 corresponding to 
-1.7 < AN,, < 2.0 (1.3 < N„ < 5.0). For an "interest- 
ing" range in £ e over which the fits are reasonably accu- 
rate, —0.1 < £ e < 0.1 is adopted. 

Before presenting our fits, it is worth reemphasizing 
that Figures 1 & 2 show clearly that D and 7 Li are 
good baryometers, sensitive to 77, while 4 He is a good 
chronometer, sensitive to 5, as well as a good leptomc- 
ter, sensitive to £ e . 



Helium-4 




FIG. 3: BBN-predicted isoabundance curves for 4 He in the 
expansion-rate parameter (S), baryon density parameter (77) 
plane. From bottom to top, Y P = 0.23, 0.24, 0.25. The solid 
curves are the BBN-predicted results, while the dotted curves 
are our fits (see the text). 
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FIG. 4: BBN-predicted isoabundance curves for 4 He (solid) 
in the lepton asymmetry parameter (£ e ), baryon density pa- 
rameter (77) plane. From top to bottom, Yp = 0.23, 0.24, 
0.25. The solid curves are the BBN-predicted results, while 
the dotted curves are our fits (see the text). 



logarithmically with 77 over the range of 77 adopted here 



(4 < 7710 



< 



Thus, while for SBBN (5 = 1 and £ e = 



0) to a very good approximation Y = Y(jj) ex In 77, we 
prefer to find a linear relation between Yp and 77. Moti- 
vated by simplicity, the following linear fit for Y versus 
77 agrees with the SBBN-predicted predicted abundance 
to within 0.0006 (< 0.25%) over our adopted range in 77. 



Yt 



FIT _ 



0.2384 + 0.0016t7 10 = 0.2384 + 77 10 /625. ( n ) 



We note that over the same range in 77 this fit also agrees 
with the BNT [16] predictions for Y P to within 0.0002 
(;$ 0.1%) or better. While the accuracy of this fit is cer- 
tainly not perfect, the difference is still smaller than the 
current errors in the observationally inferred primordial 
value of Yp . 

As with the Yp - rj relation, the following linear fits 
to the Yp - 5 and £ e relations work very well over the 
adopted parameter ranges (see Figures 1, 2, 3 & 4). 



y\ 



FIT _ 



0.2384±0.0006 + 0.0016t7i 



mo 



and 



0.16(5- l)-0.23f e . 

(12) 

AN, < 1, AY « 



Notice that for fixed 
0.013AN,,. 

It is convenient to rewrite this fit of Yp(t7, 5, £ e ) 
(eq. 12) in a form which will facilitate comparison with 
comparable fits to the vd(t7, 5, £ e ) and yu(m £>■> £e) rela- 
tions. To this end, we introduce 77H0, defined as follows, 



77ho = 625(Y P - 0.2384 ± 0.0006). 



Then, 



?7hc 



77x0 + 100(5-1) 



575£ e 



(13) 



(14) 



The 4 He abundance (mass fraction: Yp) is very in- 
sensitive to the baryon density at BBN, varying roughly 



The meaning of 7yn c is clear; ?7hc is the SBBN (5 = 1 and 
£ e = 0) value of 7710 corresponding to the adopted value 



of Yp. Once Yp is chosen, the resulting value of ?7hc 
provides a linear constraint on the combination of 77, S, 
and £ e shown in eq. 14. As may be seen in Figures 3 & 4, 
this fit works well (for 4 < 7710 < 8) for 0.23 < Y P < 0.25, 
corresponding to — 5 ^ rfn c < 7. 

Adopting the IT [32] value Y P = 0.2429 ± 0.0009 (see 
§ II), leads to 77^ = 2.81 ± 0.68 (fl^ c h 2 = 0.0103 ± 
0.0025). Note that if, indeed, the very small IT error 
estimate is realistic, then helium is a competitive bary- 
ometer to deuterium! The comparisons of this fit with 
the results from our BBN code are shown in Figures 3 & 
4. As these figures and eq. 14 show, tijjc, which is lin- 
ear in Yp, is very sensitive to S and £ e , but relatively 
insensitive to ?y. 



Deuterium 




FIG. 5: BBN-predicted isoabundance curves for D in the 
expansion-rate parameter (S), baryon density parameter (77) 
plane. From left to right, j/d = 4, 3, 2. The solid curves are 
the BBN-predicted results, while the dotted curves are our 
fits (see the text). 



Over the restricted range in 77 under study here, the 
deuterium abundance yu = 10 5 (D/H), is well described 
by a power law in r\ with j/d oc rj~ le . For our adopted 
range of 77, yu = yu{v) is we U n t by 



yf7 T ee 46.5t ?1 - 1 - 6 . 



(15) 



While the true j/d ~~ V relation is not precisely a power 
law, for 4 < 7/10 ^S 8, this fit is accurate to better than 
1%, three times smaller than the ~ 3% BBN uncertainty 
estimated by BNT [16]; this fit agrees with the BNT 
prediction to 2% or better over this range in r\. 

In analogy with rjn c defined by the 4 He abundance 
(eq. 12), it is useful to define a similar parameter, 77D, 
for D (including a 3% error estimate), 



Vb = ( 



46.5(1 ± 0.03) u/1 



/1.6 



2/D 



(16) 
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FIG. 6: BBN-predicted isoabundance curves for D in the lep- 
ton asymmetry parameter (£ e ), baryon density parameter (77) 
plane. From left to right, j/d = 4, 3, 2. The solid curves are 
the BBN-predicted results, while the dotted curves are our 
fits (see the text). 



Substituting the abundance estimate (§11) for primordial 
D into this equation leads to our adopted value for 77D = 
,+0.68 t n D h 2 = o.0221 + ? 0025 ^ 



6.06 



-0.53 



m 



])■ 



For the adopted ranges in 77, S, and £ e , a simple linear 
relation which provides a good fit to ?7d (see Figures 5 & 
6) is, 



t?d =77io-6(5'- !) + -£-. 



(17) 



As may be seen from Figures 5 & 6, this fit works quite 
well (for 4 < 7710 < 8) for 2 < t/d ^ 4, corresponding to 
5 ;$ rjo ^ 7. Given the restricted ranges for S and £ e , it 
is clear from Figures 5 & 6 and eq. 17 that D is a senstive 
baryometer (770 « 7710). 



Helium-3 



The dependence of the BBN-predicted abundance of 
3 He on the baryon density, expansion rate, and lepton 
asymmetry is similar to that of D but, 3 He is consider- 
ably less sensitive to them than is D. Indeed, over the 
parameter ranges adopted here, the 3 He abundance is 
well fit by 7/3 oc ?/p> 35 . As a result, along with its more 
complicated post-BBN evolution, 3 He provides comple- 
mentary but less compelling constraints than does D. In 
our subsequent analysis, 3 He is set aside in favor of D. 




the adopted range in 77, 



FIG. 7: BBN-predicted isoabundance curves for 7 Li in the 
expansion-rate parameter (S), baryon density parameter (77) 
plane. From left to right, j/Li = 3, 4, 5. The solid curves are 
the BBN-predicted results, while the dotted curves are our 
fits (see the text). 
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FIG. 8: BBN-predicted isoabundance curves for 7 Li in the 
lepton asymmetry parameter (£ e ), baryon density parameter 
(77) plane. From left to right, 7 Li, j/li = 3, 4, 5. The solid 
curves are the BBN-predicted results, while the dotted curves 
are our fits (see the text). 



nl FIT _ %0 



(18) 



While this fit predicts slightly smaller lithium abun- 
dances compared to those of BNT [16], the differences 
are at the 5-8% level, small compared to the BNT un- 
certainty estimates as well as those of Hata et al. [17] 
(-10-20%). 

In analogy with 77D and 7/Hc defined above, we intro- 
duce 77Li (allowing for a 10% uncertainty), defined by, 



77 Li = (8.5(1 ±0.1)y u ) 1/2 . 



(19) 



Using the Ryan et al. [34] estimate (§11) y u = 1.23to' 34 



3^ 2 3+o.4 4 



16- 



7/ Li = iS.^l Sis (^b^ = °- 0118 -o.ooio)> while the Boni- 
facio et al. [35] result of yu — 2.19^q;38 gives r/u = 
4.31 ig;^ (n^h 2 = 0.0157 ± 0.0017). 

A simple, linear relation for rju as a function of 77, 5, £ e , 
which works reasonably well over the adopted parameter 
ranges (see Figures 7 & 8) is, 



Vu = r) W ~ 3(5 - 1) 



76 



(20) 



As Figures 7 & 8 reveal, this fit works well (for 4 < 
Vw < 8) for 3 < i/Li S 5, corresponding to 5 < rju ^S 7. 
We note that this fit breaks down for yu < 2 (i%i < 4). 
In addition, Figures 7 & 8, along with eq. 20, show that, 
as is the case for deuterium, lithium can be an excellent 
baryometer (for the restricted ranges of 5 and £ e , rju ~ 

mo)- 



IV. APPLICATIONS 

In this section the application of our simple, linear fits 
(among 77, 5, £ e , and vd, YLi, Yp) is illustrated by con- 
sidering several standard and nonstandard BBN options. 

A. SBBN 



D. Lithium-7 

As with D, the 7 Li abundance 2 is well described by 
a power law in 77 over the range in baryon abundance 
explored here: yu = 10 10 (Li/H) oc ?7 2 . The following fit 
agrees with the BBN predictions to better than 3% over 



2 It is common in the astronomical literature to present the lithium 
abundance logarithmically: [Li] = 12+log(Li/H) = 2 + log(yLi). 



As a first application of our simple fits to the pre- 
dicted primordial abundances, consider the case of SBBN 
(5 = 1, £ e = 0). In this case the predicted abun- 
dances depend on only one adjustable parameter, the 
baryon abundance parameter 77. For SBBN we may 
adopt the WMAP-inspired value [8] of 7710 = 6.14 ± 0.25 
(n B h 2 = 0.0224 ± 0.009). For 5 = 1 and £ e = we 
expect, 



SBBN : 7/d = 77Hc = 7/Li = 7710. 



(21) 



It is clear from Figure 9, using the abundances identified 
in §11, that eq. 21 is not satisfied (770 rs 77 10 ^ 77 Hc ~ 
77Li). However, SBBN does work for deuterium: 77D = 
6 - 06 ia53 ~ Vw = 6.14 ± 0.25. In addition, we note that 




TllO 



FIG. 9: The baryon abundance parameters (see §111) corre- 
sponding to the D, 4 He, and 7 Li abundances adopted here 
(from §11), and from WMAP [8] 



B. Non-Standard Expansion: S ^ 1, £ e = 

The previous section has reminded us of the tension 
between the observationally-inferred primordial abun- 
dances of D, 4 He, and 7 Li; see, e.g., [38]. While it is 
not unlikely that the resolution of these conflicts may 
be found in systematic uncertainties of the astronomy 
(e.g., corrections to the derived helium abundance and/or 
lithium depletion/dilution via mixing), they may be pro- 
viding hints of new physics. Since the BBN abundance of 
He is sensitive to the early-Universe expansion rate and 
that of D probes the baryon density, a combination of 
the two abundances permits us to investigate constraints 
on models with non-standard expansion rates [11, 39]. 
Setting aside 7 Li for the moment and only using D and 
4 Hc, we may use equations 13 and 16 to find 



5-1 



?7Hc - Vd 
106 



(22) 



and, 



for the WMAP baryon abundance the SBBN-prcdictcd 
value of the 3 He primordial abundance, 1.1±0.1 x 1CP 5 , is 
in excellent agreement with that inferred from the study 
of Galactic Hn regions [19]. The problems arise for 4 He 
and 7 Li. 

For 4 He, the recent Izotov & Thuan [32] determination 
of the primordial mass fraction leads to t^hc = 2.81 ±0.68 
(§111 A), which is ~ 6c below the WMAP baryon density 
parameter. While Izotov & Thuan have attempted to 
account for some of the potential systematic errors in 
their Yp determination, it may well be that they have 
underestimated the residual error; see, e.g., Olive and 
Skillman [33]. If not, this disagreement may provide a 
hint of new physics beyond the standard model (5 5^ 1 
and/or £ e ^ 0). 

As may be seen from Figures 1, 2, 4, & 7, along with 
cqs. 16 & 19, the BBN-predictcd abundance of 7 Li is 
closely tied to that of D, yLi « 4.3(t/d/2.6) _5//4 . Forynin 
its observed range (and/or for 7710 in the WMAP range) , 
y Li w 4 - 5 ([Li] « 2.6 - 2.7). For the WMAP baryon 
abundance the SBBN prediction is y Li = 4.44t°'57, to 
be compared with the halo star estimate from Ryan et 
al. [34] of yLi = 1.23^oi6> or with the globular cluster 
estimate from Bonifacio et al. [35] of yLi = 2.19^Q 3 g. 
The SBBN-predicted 7 Li abundance is much higher, by 
factors of ~ 2 — 4, than those inferred from the stud- 
ies of the oldest halo and/or globular cluster stars in 
the Galaxy. It will be seen in more detail below that 
it is unlikely that any of the "simple" nonstandard ef- 
fects (S 7^ 1, £ e 7^ 0) can resolve this conflict (since the 
7 Li and D abundances are so tightly coupled). It may 
well be that for these oldest stars in the Galaxy, mixing 
surface material with the lithium-depleted interior has 
reduced their observed, surface lithium abundances from 
the primordial value [36, 37]. 



Vw 



lOOr/p + 6?7hc 
106 ' 



(23) 



It is clear from eq. 23 that D is the dominant baryometer. 
While it may appear from eq. 22 that D and 4 He make 
comparable contributions to constraining the expansion 
rate, recall that ?7ho is linear in Yp while 770 varies only 
as a fractional power of vd- Thus, whereas reconciling 
Yp with yn requires an increase in Yp of AYp w 0.006 
(~ 2.5%), reconciling D with 4 He would require that vd 
increase by a factor of ~ 4 — 5 (and, the resulting tj-q 
would then be in conflict with the WMAP determination 
of 7710). From eq. 22, for AN„ -C 1 and with vd (??d) 
fixed, AY « 0.014AN„. Note that depending on the 
quantity held fixed (7710 or tjd), the coefficients in the 
AY - AN„ relations differ slightly. 

For the values adopted in §IIIA,B (and their uncer- 
tainties), 77 D = 6.06±g;|| and r] Hc = 2.81 ± 0.68, we find 



mo 



-0.64 
-0.50 



S = 0.969: 



0.009- 



(24) 



If the non-standard expansion rate factor is associ- 
ated with an equivalent number of neutrinos, AN„ = 
-0.37to:!i- These results, along with the D and 4 Hc 
isoabundance curves, are shown in Figure 10. Notice that 
although AiV„ is now closer to zero than in, e.g., [9], the 
smaller adopted uncertainty in Yp results in a larger, 
^3 — 4(7 discrepancy. However, these BBN results for 
the baryon density and the expansion rate are in excellent 
agreement with those from WMAP which, while sensitive 
to the baryon density parameter is relatively insensitive 
to the expansion rate factor; see [9]. 

Although the above combination of 7710 and S can rec- 
oncile the inferred primordial abundances of D and 4 Hc, 
it cannot resolve the lithium problem. 



Vu 



10377D + 377Hc 

106 



5.97; 



-0.66 

-0.52- 



(25) 
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FIG. 10: BBN-predicted isoabundance curves for D (dashed) 
and 4 He (solid) in the expansion-rate parameter (S), baryon 
density parameter (r/) plane. The point with error bars cor- 
responds to the adopted abundances for D and 4 He (£ e = 0) ; 
see the text, eq. 23. 



Accounting for the theoretical uncertainty in the BBN 
lithium production and in our fit, 



>'Li 



= 4.20 



+ 1-07 . 
-0.81 ' 



[Li] P = 2.62±° 



0.10 
09- 



(26) 



This conflict is shown in Figure 11 where the D and 4 He 
constraints on 7710 and S are shown along with isoabun- 
dance curves for 7 Li. It is clear that yLi < 3 is strongly 
disfavored. 



C. Non-Zero Lepton Number: £ e 7^ 0, S — 1 

In an attempt to relieve the tension between the BBN 
predicted and observed deuterium and helium-4 abun- 
dances, a non-zero lepton number (v e — D e asymmetry; 
£ e 7^ 0) is an alternative to the non-standard expansion 
rate explored in the previous section. An excess of v e 
over v e will drive down the pre-BBN neutron-to-proton 
ratio, leaving fewer neutrons to be incorporated in 4 He. 
The resulting, lower, BBN-predicted 4 He abundance will 
be closer to that observed. At the same time, the effect of 
an e-neutrino asymmetry on the abundances of D and/or 
7 Li is small. 

As for the case of S ^ 1, we set aside 7 Li and use D 
and He to constrain 7710 and £ e - For S = 1, 



>/io = 



115t7e 



Va c 



116 



(27) 



and, 




FIG. 11: BBN-predicted isoabundance curves 7 Li in the 
expansion-rate parameter (S), baryon density parameter (77) 
plane. The numbers next to the curves are for yLi- The point 
with error bars corresponds to the adopted abundances for D 
and 4 He (£ e = 0) as in Figure 9. 



L = 



Vu - Vn c 
145 



Using the D and 4 He abundances adopted here, 



Vw = 6.04l«j 7 



t - 022+ 0007 

Se — u - uzz -0.006 



(28) 



(29) 



These results, along with the D and 4 He isoabundance 
curves, are shown in Figure 12. 

A "small" lepton asymmetry, £ e sw 0.02, which, how- 
ever, is very large compared to the baryon asymmetry 
(~ 77 < 10~ 9 ), can reconcile the observationally inferred 
primordial D and 4 He abundances with the predictions 
of BBN. But, this lepton asymmetry will not resolve 
the conflict between the BBN-predicted and observed 7 Li 
abundances since 



'/Li 



142T7D + 37/He 



145 



and, 



yLi 



- 4 23+ 108 



[Li] P 



5 qq+ - 67 

,J - yt '-0.52 



■ 2 6T f01 ° 

- i-UO_Q Q 9 . 



(30) 



(31) 



As is the case for a non-standard expansion rate, the 
BBN-predicted abundance of 7 Li is very tightly tied to 
that of D, and a non-zero lepton number cannot resolve 
the conflict between theory and data. This is illustrated 
in Figure 13, where it is clear that here, too, yLi < 3 is 
strongly disfavored. 
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FIG. 12: BBN-predicted isoabundance curves for D (dashed) 
and 4 He (solid) in the lepton asymmetry parameter (£ e ), 
baryon density parameter (rf) plane. The point with error 
bars corresponds to the adopted abundances for D and 4 He 
(S = 1); see the text, eq. 28. 




FIG. 13: BBN-predicted isoabundance curves 7 Li in the lep- 
ton asymmetry parameter (£ e ), baryon density parameter (if) 
plane. The numbers next to the curves are for yLi- The point 
with error bars corresponds to the adopted abundances for D 
and 4 He (5* = 1) as in Figure 11. 



D. Non-Zero Lepton Number (£ e / 0) And 
Non-Standard Expansion Rate (S / 1; N„ / 3) 

In §IVB it was seen that, in the absence of a non- 
zero lepton number, the D and 4 He abundances severely 
restrict deviations of the early universe expansion rate 
from its standard value. When associated with the effec- 
tive number of neutrinos, this constraint (for the choice 
of D and 4 He abundances adopted here) is nearly 4a away 
from the standard value of N„ = 3. If the LSND result 
[40] is interpreted in terms of a 4th, "sterile" , neutrino, 
the mixing of such a neutrino with the active neutrinos 
would bring it into equilibrium in the early universe prior 
to BBN [41], resulting in N„ = 4 at BBN. Clearly (see 
§IIIB), the current best estimates of the primordial abun- 
dances cannot tolerate the corresponding speed-up in the 
expansion rate at BBN. However, in §IIIC it was shown 
that N y = 3 is perfectly acceptable provided that there 
is a small but significant v e — v e asymmetry. It is, there- 
fore, not unexpected that for £ e ^ 0, the BBN constraints 
on 5 (AN„) are considerably expanded (see, e.g., Barger 
et al. [9] and references therein). In implementing for 
BBN 5 7^ 1 as well as £ e ^ 0, there is a problem. Be- 
cause of concerns about the primordial 7 Li (and 3 He) 
abundances as inferred from current observational data, 
the analysis here has been restricted to employing only 
the D and 4 He abundances. But, with three free param- 
eters (?yio, 5 (or AN„), and £ e ) and only two constraints 
(r)D and ?7hc), non-BBN data (restricting 7710) is required 
to proceed further. Retaining the constraints from D and 
4 He, 5 and £ e are then functions of the baryon density 
parameter. 



590(5 - 1) = 116r7io - (H5r?D + ??Hc) 



(32) 



and, 



145£ e = 106(5 -1) + tid-tihc- (33) 

For the abundances adopted here (§11), this leads to 



590(5-1) = II677 



10 



700; 



-78 
-61' 



and 



£ e = 0.731(5-1) + 0.0224 



-0.0066 
-0.0059- 



(34) 



(35) 



This £ e - 5 relation is shown in Figure 14. 

A constraint on the baryon density is possible utilizing 
CBR data provided that the corresponding values of 5— 1 
and £ e are sufficiently small, so that the expansion rate at 
recombination, some 400 kyr after BBN, doesn't deviate 
significantly from the standard value. For the WMAP 
recommended baryon abundance [8], 7710 = 6.14, 5 ~ 
1.0203, corresponding to AN„ « 0.25 and £ e « 0.037. 
This point is shown in Figure 14. 

It is now possible to address the question of whether 
or not BBN permits N„ = 4 (5 = 1.0783). This choice 
for 5 requires £ e = 0.080 and 7710 = 6.43, only ~ la away 
from the WMAP value. This point is shown in Figure 
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FIG. 14: BBN-predicted £ e — S relation corresponding to the 
adopted D and 4 He abundances. The solid curve is for the 
central value and the dashed curves represent the la uncer- 
tainties. The point indicated corresponds to the WMAP value 
for the baryon density; see the text. 



15 where the £ e — AN„ relation is displayed. While BBN 
constraints on more extreme cases (e.g., AN„ > 2) can 
be explored (see, e.g., Barger et al. [10]), our simple 
fits begin to become inaccurate for £ e > 0.1, limiting our 
analysis here to AN„ ^ 1.2 (see Figure 15). 

For the general case studied here, the BBN-predicted 
lithium abundance depends on the baryon density (as 
well as on the D and 4 He abundances). 



rju = 0.1597/m + 0.822t7 D + 0.019?7 H c 



(36) 



Notice that here, too, the 7 Li abundance is tightly tied 
to that of D. Indeed, 



Vu = Vb+ 3[(5 - 1) - £ e ] « T7D- 



(37) 



As before, there are lithium "problems" for the two cases 
studied above. For the WMAP combination of parame- 
ters (see eq. 19), r/u — 6.01, so that yu — 4.25 ([Li] = 
2.63). For N„ = 4, r/u — 6.06, corresponding to yu = 
4.32 ([Li] = 2.64). At least for these combinations of non- 
standard physics (as well as for SBBN), there is an ap- 
parent conflict between the BBN predicted and observed 
primordial abundances of lithium. 




FIG. 15: BBN-predicted £ e — AN„ relation corresponding to 
the adopted D and 4 He abundances. The solid curve is for 
the central value and the dashed curves represent the la un- 
certainties. The point indicated corresponds to N„ = 4. 



SUMMARY AND CONCLUSIONS 



In the toolbox of a particle phenomenologist, BBN is 
a venerable and valuable tool, providing complementary 
constraints on - or hints of - new physics beyond the 
standard model. As the models change and the data are 
updated, not every phenomenologist on the street may 
have easy access to a BBN code to enable quick revision 
of previous bounds or investigation of new hints. It may, 
then, be of value to have simple fits to the BBN-predicted 
abundances as functions of the key variables (baryon den- 
sity, expansion rate or neutrino number, neutral lepton 
asymmetry). For interesting but limited ranges of the 
parameters such fits have been presented here and ap- 
plied to several simple examples. We have seen that the 
observed D and 4 He abundances, while apparently in- 
consistent with SBBN (modulo systematic errors in the 
data) , can be reconciled if non-standard expansion rates 
and/or non-zero lepton number are allowed. It was also 
shown how non-zero lepton number allows one - or more 
- sterile neutrinos while maintaining the consistency of 
BBN. Finally, by revealing the clear connections between 
the BBN-predicted abundances of D and 7 Li, we showed 
that none of these non-standard physics solutions can rec- 
oncile the observed and predicted lithium abundances. 
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